Abstract 

A possibility of geophysical measurements using the large scale laser in- 
terferometrical gravitational wave antenna is discussed. An interferometer 
with suspended mirrors can be used as a gradiometer measuring variations 
of an angle between gravity force vectors acting on the spatially separated 
suspensions. We analyze restrictions imposed by the atmospheric noises 
on feasibility of such measurements. Two models of the atmosphere are 
invoked: a quiet atmosphere with a hydrostatic coupling of pressure and 
density and a dynamic model of moving region of the density anomaly 
(cyclone). Both models lead to similar conclusions up to numerical fac- 
tors. Besides the hydrostatic approximation, we use a model of turbulent 
atmosphere with the pressure fluctuation spectrum ~ f~ 7 ^ 3 to explore the 
Newtonian noise in a higher frequency domain (up to 10 Hz) predicting the 
gravitational noise background for modern gravitational wave detectors. 
Our estimates show that this could pose a serious problem for realiza- 
tion of such projects. Finally, angular fluctuations of spatially separated 
pendula are investigated via computer simulation for some realistic atmo- 
spheric data giving the level estimate ~ 10 -11 rad-Hz _1,/2 at frequency 
~ 10 -4 Hz. This looks promising for the possibility of the measurement 
of weak gravity effects such as Earth inner core oscillations. 

Introduction 

A possibility of geophysical measurements with large scale laser intcrfcromctrical 
gravitational wave antenna [El || was discussed in |3j [| . It was supposed that at 
very low frequencies 10~ 3 — 10 Hz the interferometer with suspended mirrors 
can be considered as an angular gravity gradiometer measuring variations of the 
mutual angle between gravity force vectors or plumb lines of spatially separated 
suspensions. Geophysical information could be read out as error signal of the 
feedback circuits which preserves the operational angle position of the mirrors. 

As an example of geophysical phenomenon for the measurement a very weak 
effect of the Earth inner core oscillations (one of "hot points" of modern geo- 
physics ||) was considered in 

A number of principal instrumental problems were analyzed in ^ Q (such 
as a necessity of suspension for the laser injection bench, the problem of tilt and 
shift differentiation for spherical mirrors, and others); the main part of instru- 
mental noises at low frequencies was estimated in where a positive conclusion 
was drawn in favor of a feasibility of geophysical measurements. However, envi- 
ronmental geophysical fluctuations produce a predominant noise background 
which can be a problem while realizing this program. It is clear that, for 
any gravitational device, the fundamental and unavoidable source of noises is 
the gravitational attraction of stochastically moving surrounding masses, the 
so-called "Newtonian gravity noise." Such noise produced by acoustic waves, 
propagating under the interferometer base ground, was calculated in 8 at high 
frequency / >10 Hz typical for gravitational detectors. 

In this paper, we estimate, at least at the order of magnitude, the Newtonian 
fluctuations of the mutual angle between two separated plumb lines (suspended 
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mirrors) produced by atmospheric perturbations; preliminary results were pre- 
sented in B. 

Stochastic density variations and transportation of large air masses produces 
a contribution into the gravity force on the Earth surface. Almost all processes 
of such kind are connected with the corresponding pressure changes. Thus 
they can be controlled by means of pressure monitoring. In the gravimetry 
theory, it is well known that the influence of atmospheric pressure on the gravity 
consists of direct attraction of the atmospheric mass and crust deformation due 
to atmospheric loading [ [To| . 

The both effects have the same order of magnitude but the loading effect, as 
a rule, is at least two-five times less then the direct attraction. For this reason, 
below we consider only the direct attraction effect as one of main contributions 
to the atmospheric Newtonian noises. 

As a control reference point in our analysis, we bear in mind an acceptable 
level of the noise spectrum density which allows one to register the inner core 
oscillations with amplitudes of the order of 1 m — the expected signal can 
achieve 10 -13 — 10~ 14 rad for the interferometer base equal to 3 km j|, ||. 
Suppose that the period and quality factor of the inner core of the fundamental 
mode are r = 3.3 h and Q = 10 2 — 10 3 , respectively Then one can take 
the measurement time r m = 10 6 sec which results in the following angular noise 
spectrum limit Sctf = 10~ n rad-Hz -1 / 2 . 



1 Plane atmosphere 

A very rough and simple Earth atmospheric model is a half-space filled by 
air stratum of the homogeneous density p of the altitude h. A change in the 
gravity acceleration Ag caused by perturbation of the atmospheric density can 



be estimated as follows |12 



. 2nG . . . 

Ag z = — Ap, (1) 

where Ap is the atmospheric pressure variation, g is the gravity acceleration, 
G is Newtonian gravity constant and z-axis is directed upward. The factor of 
admittance between the gravity and pressure k p — 2irG/g = 0, 427/^Gal/hPa is 
widely used in geophysics and is in a fair correspondence with experimental data. 
One can easily get (|]) starting with the surface field of the stratum g z = 2irGph 
and using a hydrostatic coupling between the pressure and density p — pgh if 
one supposes very slow variations. 

Although this formula has an integral sense, one can likely use it in a spectral 
form by substituting the corresponding spectral densities of the gravity and 
pressure variables, having in mind the fact that the changes of pressure are 
mostly quasistatic processes, i.e. the factor of admittance does not depend on 
frequency. 

After this, the relative angular deflection of two plumb lines located at the 
ends of the interferometer's base can be evaluated phenomenologically, in view 
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of the following expression: 



<fo>^^f = £&f b(/)K/A/, (2) 

where L is the arm's length, L p is the scale of the pressure anomaly domain 
(cyclone size, etc.), p(f) is the spectral density of the presure amplitude in the 
bandwidth Af and £ <C 1 is a coupling factor between vertical g z and horizontal 
g x components of gravity variations. 

To evaluate the £ factor, we suppose that pressure changes along the x- 
direction so slowly, that ([!]) remains to be valid and it leads to the condition 
h <C L p . Then using the potential character of the gravity field, i.e., rot g = 0, 
one can take dg x /dz = dg z /dx, or Ag x /h ~ Ag z /L p . This results in the 
estimation 

Ag x h 
Ag z L p 

which might be used while calculating in view of formula (^|). 

The pressure time variation background was investigated by many meteoro- 
logical observatories in different places. In spite that it is a local characteristic, 
there are typical common features in the spectrum density such as the diurnal 
(1.17-10~ 5 Hz) and semidiurnal (2.33-10 -5 Hz) peaks, a general upward slope at 
very low frequencies, etc. For numerical estimations at the order of magnitude, 
we used data collected by the Mizusawa Astrogeodynamic Observatory (Japan) 
from 1986 to 1987. 

Available data jl3| can be approximatelyly separated into three following 
spectral domains with the dominant spectral amplitudes, namely, 

(a) (2 • 1CT 5 - 2.5 • 1CT 5 ), \p(f)\ = 9 • 10 2 mBar/Hz 1/2 

(b) (1,5- 1CT 5 - 8- 1CT 6 ), \p(f)\ = 1,9 • 10 3 mBar/Hz 1/2 

(c) (3-1CT 6 - 5-1CT 7 ), \p(f)\ = 9,5 • 10 3 mBar/Hz 1/2 . 

This demonstrates that a transition from time-scale of several hours (a) to 
time-scale of several days (b) is accompanied by the increase in the pressure 
fluctuations, in average, of one order of magnitude. One can see that the domain 
of several hours (a) overlaps with the tidal semidiurnal spectrum peak, which 
can be removed while data processing, so it is reasonable to take for estimation a 
decreased (but still dominant) extrapolation value 1-10 2 mbr/Hz 1 ^ 2 for the times 
3-1-4 hours (inner core polar mode). Then taking as typical average pressure 
anomaly the (cyclone) size L p = 500 km and the effective atmosphere altitude 
h=10 km, one obtains the following result for the standard angle deviation 
between arm's mirrors 



<Sa>~ ^ -2 \p(f)\VAf ~ 5 • 10- 12 VA/ rad. (3) 
9 L p 



This estimation does not exceed the critical noise density mentioned in Intro- 
duction although it shows difficulties in the detection of the inner core motion 
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because the angular noise is at the border of the expected gravity effect. How- 
ever, it is supposed that a correction for the gravity atmospheric noise could be 
carried out if one would precisely control pressure variations in the location of 
front and end mirrors. 

2 Cyclone dynamic effect 

An obvious lack of the previous consideration is an uncertainty of the £ factor 
which cannot be regorously found within the framework of the homogeneous 
static model of the atmosphere. To avoid this uncertainty, it would be useful to 
consider a more complex dynamic model where some spatially limited domain 
inside the atmospheric half space, called below as "cyclone" (characterized by 
perturbed parameters of density and pressure) moves along the interferometer 
base. Then a direct calculation of its Newtonian attraction permits one to get 
a more reliable estimation of mutual angular deflections of the interferometer 
arm's mirrors. 

It is convenient to take the cyclone form as an oblate ellipsoid of rotation 
because its gravitational field can be expressed in terms of elementary functions 
[ p"4[ . After this, the problem is formulated as follows. The cyclone (ellipsoid 
of rotation) with its plane of symmetry on the Earth surface moves along the 
x axis where the interferometer arm is located (figure 1). Only upper half of 
the ellipsoid has a physical sense representing the cyclone, so that an effective 
gravitational field along the x- axis must decrease twice. A vertical size of the 
cyclone (a minor ellipsoid semiaxis c) is approximately equal to the atmosphere 
altitude in z direction (c < h), however, it is much less then a horizontal size of 
the cyclone, a major semiaxis a (a c ~ h). Let us take a Cartesian coordinate 
system in the cyclone center in such a way that the current value x is just the 
distance between the cyclone and interferometer centers. The positions of the 
front and end suspended mirrors are (x—L/2) and (x+L/2), with L representing 
the interferometer arm length. The cyclone pass through the interferometer base 
corresponds to the x variation in the limits (+00; —00). 

For a homogeneous ellipsoid, its gravity field (acceleration) along the x axis 
grows linearly with the distance from the center inside the figure (x < a), while 
outside (x > a) it can be described by the following expression 



g x = -2tt Gp—p\r) 1 \/i] 2 - 1 - 77 arcsinr? 1 j , (4) 

where A 2 — a 2 — c 2 , t] = x/A and p is an excess (or deficit) of air density over 
the average density of unperturbed atmosphere. 

In principle, formula (^), in combination with the corresponding formula 
inside the figure, gives a possibility to determine dynamic perturbations of mir- 
ror's plumb lines when the cyclone is penetrating through the base. However, 
it is much more convenient to operate with the field of a thin ellipsoid shell to 
avoid a necessity of the transition between two different field expressions. 
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To get formula for the shell field, it is enough to accept c = ka and then to 
differentiate g x (|4|) as a function of the variable a. This procedure results in the 
following formula: 

dg x =AnG P k — = da (5) 
Xyjx 2 — a 2 (l — k 2 ) 

Expression (^) presents the gravitational field outside the shell (x > a) , mean- 
while, the internal field is equal to zero, the factor da is the shell's thickness in 
the cross point with the x axis. Then the field of inhomogeneous cloud at the 
points of x axis can be calculated in terms of the integral 

g x = 4*Gk [ X , (6) 

Jo X\Jx 2 — a 2 (l — k 2 ) 

where the function p(a) defines the cyclone density distribution (the maximum 
horizontal size of the cyclone, hereinafter, noted as ao). 

The mutual angular deflection of mirror's suspensions located at the ends of 
the interferometer base is given by the simple combination of the function g x {x) 
in two particular points of the x axis: 

g x {x + L/2)-g x {x-L/2) 
5a = . (7) 

Factor 2 in the denominator reflects the effect of the upper part of the ellipsoid 
cloud only, i.e. the domain where z > 0. For transition to the time domain, one 
has to substitute x = vt, where v is the cyclone velocity. 

To illustrate a general type of the angular deflection dynamics, formulas (||) 
and (|7]) were calculated numerically for several artificially selected laws of the 
density variation over the cyclone. Four types of the density distributions were 
considered, namely, "step," "circle," "harmonic" and "Gaussian" (see figure 2, 
curves la-Id) with the same normalization of mean value and variance: 

P ( X )kx> dx =(./2)- p alc, f P (x)kxUx = (,/2)- p at, 

Corresponding plots are presented in figure 2. Density distributions used for 
calculations are presented by curves la-Id and curves 2a-2d show the dynamics 
of relative variations of the mirror's plumb lines. For numerical calculations, it 
was used v=l m/sec and the other parameters were ao=300 km, c=10 km, 
p=0.02 kg/m 3 and L=3 km. 

Figure 2 demonstrates a strong dependence of the mutual angular deflection 
on the density distribution inside the cyclone. Common features are peaks 
of deflection at the cyclone base borders and a suppression of perturbation 
inside the cyclone. The amplitudes of peaks depend on the density structure 
of the cyclone cloud. The homogeneous cloud produces the biggest jump of 
deflection (case a). The Gaussian distribution minimizes the jump effect (case 
d). However, the best compensation of perturbations inside the cyclone takes 
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place just for the homogeneous model (here the mutual deflection is more close 
to zero then in the other cases). 

For the homogeneous model (case a), the integral in (^) can be calculated 
analytically. Results of such calculation within the first order accuracy of the 
parameter (c/ao) <C 1 read 

g x (x)~—ir 2 pG — x for \x\ < ao; 
a 

(dg/dx) ~ ir 2 pG — for |x| < ao: 
a 

Ag = g x (x + L/2) - g x (x - L/2) ~ g x \x=a —• 

a 

For the mutual angular deflection, one has 

2g 2g a 4 p a W 

The adiabatic coupling between the pressure and density variations (p/po) — 
(Ap/po) was used to obtain relationship (||), where v 2 ~ Po/po is the sound 
velocity in the normal atmosphere and Ap is the pressure jump outside the 
cyclone. By comparing (||) with phenomenological formula one can see that 
the formulas are similar. For c ~ h ~10 km, (h/v 2 ) — l/g (v s —33Q m/sec) one 
obtains £ ~ 1 for the rough homogeneous model, but it decreases up to 0.01 for 
more smooth cyclone shapes. 

From figure 2 it is clear that a more favorable cyclone structure can be pre- 
sented by a cloud with Gaussian fronts and long plane central zone. Then jump 
perturbations corresponding to the cyclone front passing through the interfer- 
ometer base are suppressed as well as perturbations in the central part of the 
cloud. In general, large sharp perturbations (corresponding to relatively rare 
events) can be obviously controlled by barometric measurements. At quiet at- 
mospheric conditions, the usual daily variation is of the order of (1 -j- 10) mBr, 
that keeps the effect of Newtonian atmospheric attraction, in view of (||) , at the 
level of 10~ n ~ 10~ 12 rad. 

Nevertheless, it is by two-three order of magnitude larger that the inner core 
oscillation effect (~ 10~ 13 -r 10~ 14 rad). However, (g) gives the integral value 
of the atmospheric gravity effect and one has to evaluate its spectrum noise 
density in the domain of interest at / ~ (10~ 4 -^10~ 3 ) Hz. 

Let us consider Poissonian flux of cyclones with average number n per time 
unit. Excluding the cases of recovering, one can take n < v/2a. The well-known 
formula for the power spectrum of the Poissonian flux defined by the spectrum 
of individual events reads 

SM) = n\Sa(.f)\ 2 , (9) 

where 6a(f) is the Fourier component of the angle deflection. It is clear a priori 
that the spectrum maximum takes place near / *~ v/a = 10~ 6 Hz independently 
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of the cyclone density distribution. The spectrum behavior at high frequencies 
strongly depends on the cyclone density structure. The main role is played by 
the peak amplitude at the border x = vt = a. For the thin ellipsoid shell, 
the peak amplitude is AirGpda and it gives the elementary gravity variation 
spectrum 

ApGv 2naf 

dgx(J) = — 7— cos . 

/ v 

It has to be integrated over the cyclone cloud in order to get the total gravity 
spectrum 

*c/)-jf*.(/)-^,(¥), m 

i.e. it is essentially determined by the cloud density spectrum p. To obtain 
the mutual angular spectrum of plumb lines, one has to multiply (|ic| ) by the 
transfer function of the differential link (Q). Numerical calculations for the four 
selected cyclone clouds are presented in figure 2 by curves 3a-3d. 

Let us consider the homogeneous cyclone (case a) in more detail. The in- 
tegrand in ( |Io| ) is calculated analytically and the final result for the angular 
deflection power spectrum reads 

Mfl-*2p£ *.(?=/) (") 

At / > v/L ~ 10~ 4 Hz, the spectrum oscillates and falls down as l// 4 and 
at low frequencies S a ~ l// 2 - At frequency / ~ 1CP 4 Hz (interesting for us), 
the angular spectrum density equals to 2 • 10~ 9 rad-Hz -1 ^ 2 , that is two order of 
magnitude larger than the acceptable level. However, it follows from the very 
sharp density jump in the homogenous cyclone model. For more smooth cyclone 
structure, the angular spectrum at high frequencies goes down much faster: 

S a ~ 1/f 3 (3b); 

S a ~ 1/f 4 (3c); 

S a ~exp-/ 2 (3d). 

From figure 2 one can see, that at frequency 10 -4 Hz the amplitude of 
angular noise is already less than the critical level 10~ n rad-Hz -1 ^ 2 (3c) and 
then it becomes negligible (3d). 

Analysis presented in sections 1 and 2 elucidates particular characteristics of 
the free mass interferometer such as a gravity gradiometer one. Being a gravity 
device, it is affected by all environmental movements of masses. But its transfer 
function suppresses the influence of large scale homogeneous gravity sources by 
the factor of the ratio of the base length to the source scale. The most effective 
sources would be those which generate the gravity field with a space scale equal 
to the interferometer base (or less) . Such sources might have a typical variation 
time of the order of tens seconds for v > 10 m/sec, so it would be high frequency 
sources with respect to the main geophysical processes with times equal to hours 
and longer. It is clear also that, by suppressing the Newtonian environmental 
noise, the gravity gradiometer suppresses equally low frequency (quasistatic) 



7 



signals arising from the objects (including the inner core oscillations). However, 
it cannot be a serious objection if instrumental (intrinsic) noises of the device are 
enough small to permit measurements of such very weak effects. It seems that 
modern advanced technology of gravitational interferometers under construction 
could meet these requirements. 



3 Turbulent flux efect 

The previous consideration was carried out at the condition of the hydrostatic 
atmosphere or, at least, slow laminar air currents. However, it is known that in 
upper troposphere strata the Reynolds number can be very high Re> 10 6 and 
appearance of turbulent currents is a typical event. The rigorous calculation 
of the turbulent atmosphere gravity effect is complicated due to an uncertainty 
in constructing the models adequate to our problem, many unknown param- 
eters and general complexity of the theory of turbulent atmosphere. For the 
estimation of the order of magnitude of the gravity noise produced by the air 
turbulence, we hope that a simple phenomcnological consideration based on the 
Kolmogorov-Obukhov (K-O) law |15| can be useful. This law fixes relative vari- 
ation of velocities of vortex zones with different scales in the regime of developed 
turbulence. 

Let a be some maximum size of the turbulence zone with an average ampli- 
tude of the velocity variations V. Then according to the K-0 law, a velocity 
standard v\ at the small scale A < a obeys to the formula 

A XV3 



vx^Vl-) . (12) 



For each A, one can define the frequency of the velocity fluctuation as / = V/ A. 
Then ( |l2[ ) can be rewritten as 

v x = V i / 3 (fa)- 1 / 3 . (13) 



Following [|16| the pressure standard Ap is introduced in (10) through the hy- 
drodynamic law at the scale A: 

while the variance < (Ap) 2 > is composed by contributions of all turbulent 
scales A, or corresponding frequencies /(A) from a minimal scale, say, Ao up to 
the current value A, i.e. 

<(Ap) 2 >= S p (f)df -const- / S p (f)df. (14) 

Jf(\) Jo 

Here S p (f) is the spectral power density of the pressure variations and the 
approximation Aq — 0, /(Ao) — oo was used. 
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Differentiating (|14|) over frequency and taking into account ( |13[ ) , one obtains 
the well-known spectrum of developed turbulence 

S p (f)^^) 2 V^a-^f-y\ (15) 

The frequency domain where the spectrum (^) is valid reads 

- = / 2 </</ 3 = -Re 3/4 . 
a a 

At frequencies / < f% (the domain of slow hydrodynamic changes), the pressure 
variations conventionally obey to the AP ~ / _1 law (velocity fluctuations v\ 
are considered as independent), so the pressure spectrum in this domain reads 

S p (f)c^fl=S f- 3 . (16) 

Here the constant So can be determined by the condition of uninterrupted spec- 
tra in the point fi for both approaches under consideration. It is obvious also 
that the lowest conceivable frequency associated with the limited turbulence 
zone is f\ = V/ R ei where R e is the Earth radius. Below this point, the spectrum 
of pressure variations remains to be uncertain at least within the framework of 
the approach used. 

This simple model of the turbulent atmosphere pressure spectrum does not 
take into account all complex processes which could take place. So, experimental 
observations of thunder storms in the frequency domain (10~ 2 -j- 1CP 4 ) Hz gives 
the spectrum indices distributed, in average, between 2 and 4.25 and centered 
in the point 3.5 jlTfl , i.e. larger than the K-0 index 7/3. On the contrary, at 
frequencies below 10~ 4 , experimental data better agree with index 5/3 ]l8|j , i.e. 
low frequency fluctuations grow slowly than it is predicted by the flicker type 
noise with the spectrum index 3, and so on. 

Nevertheless, this model at the quiet atmosphere conditions corresponds 
to the real experimental data by the order of magnitude. In particular, it was 
successfully used for calculating the atmospheric load deformations measured by 
strain meters at the Obninsk Seismic Station (see [^6|). In figure 3 (adopted from 
[fL6|), one can compare the model spectrum calculated for the turbulent zone 
with kinetic parameters a = L = 18 km, V = AU = 3 m/sec and average density 
Po = 0.55 kg/m 3 (corresponding altitude h = 8 km) with the experimental 
pressure spectrum measured at the Obninsk Seismic Station. One can see a 
satisfactory coincidence, at least, by the order of magnitude. 

Now we discuss the gravity noise associated with a turbulent atmospheric 
process. For this, let us consider the plane atmosphere with effective altitude h 
filled by the turbulent flux as it was considered in section 2. The estimation of 
gravity variations on the Earth surface is taken as follows 

Ag ~ 2-KGhAp. 
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The average density changes Ap over the total scale of turbulent zone can be 
replaced by the average pressure standard deviation using the adiabatic law, i.e. 

Ag ~ 2nG^A P . (17) 
P 

The mutual plumb lines deflection has to be a function of the pressure gradient, 

*--h ? £ L -h™f{&)* ,18) 

(In this approximation, we suppose that the standard pressure deviations do not 
depend on x). To evaluate the pressure gradient along the interferometer base, 
one can employ again the K-0 law or use direct experimental data provided by 
some spatial barometric grid. Below we consider both methods. 

One can express the pressure gradient along a turbulence cell (A ~ L) in 
terms of the spatial derivative of the local relation p = p (v\/2) 

dp dv x T v x dv x 

d~ x = p » Vx ^ = Lpo T^- (19) 

If v\ is not known, one can take it from the K-0 law, v\ = V(Lfa) 1 / 3 , and then, 
in view of the linear dependence of the velocity on x, i.e. (v\/L) ~ (dv\/dx), 
one obtains 

l^»(l) 2 (!) 2 ' 3 ~'»- 2 ^/ tm , (20) 

where the parameters in figure 3 used for numerical estimation are a=18 km, 
V=3 m/sec and the interferometer base L=3 km. 

This very approximate estimation is supported by experimental data pro- 
vided by the barometric space pattern covered several hundred kilometers around 
the Brussels Jl0| . Measured value of the pressure gradient along the East-West 
direction was in average 0.037 hPa/km (and 0.085 hPa/km in the Nord-South 
direction) . This (by the order of magnitude) is in a satisfactory agreement with 
theoretical calculation. Below, we use the intermediate value 5 • 10 -2 hPa/km 
for estimations. 

To transform formula (|l8| ) for the relative angular variations of arm mirrors, 
one can substitute (po/po) = v$ and use the K-0 spectrum jl^ ) for estimating 
the standard deviation Ap. It results in the formula for angle deflection 

< Sa(f) > ~ ——^■^-L < S p (f) > 1 / 2 y/Zf. (21) 
pgvj ox 

For numerical estimations, let us take the following parameters: v s =330 m/sec, 
the average pressure magnitude p = 10 3 hPa, the pressure jump at the arm ends 

Ap L = ^ L = 0.05 hPa/km ■ 3 km = 0.15 hPa, 
ox 
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and the pressure spectrum density 

S p (f = 0.001Hz) ~ 10 3 Pa 2 /Hz. 

As a result, one obtains < a(f) >~ 2 • 10~ 14 rad-Hz~ 1/2 . The noise is proved 
to be small enough to be neglected in measurements of the Earth core gravity 
oscillations. 

It is worthy to estimate a residual Newtonian angular noise in the fre- 
quency domain typical for gravitational wave detectors, i.e. above 10 Hz. 
Following the K-0 law the numerical result must be multiplied by the fac- 
tor (10 Hz/0.001 Hz)" 7 / 5 = 2.5 • 10~ 6 , that yields < a(10 Hz) > ~ 5 • 
10~ 20 rad-Hz -1 / 2 (the frequency 10 Hz just corresponds to the upper boundary 
of validity of the K-0 approach if Rc> 10 6 , V=3 m/scc, a=20 km). 

For VIRGO superattenuator, a length of suspension is of the order of 10 
m, it provides the interferometer optical arm variations on the level of 5 • 
10~ 17 cm/Hz 1 / 2 . This means that there exists a deformation noise at 1.7 • 
10~ 22 Hz -1 / 2 , which might create some problems for advanced VIRGO project. 



4 Simulation of plumb lines deflection 

To test a likelihood of our estimations, we carried out a simple computer simu- 
lation of the angular dynamics of two separated plumb lines (pendula) in a vari- 
able pressure field. Real experimental data were kindly provided by Japanese 
Weather Association. These data were collected by the Meteorological National 
Geographical Institute at the area close to Tsukubo Scientific Center. The area 
20x20 km was covered by 10 meteostations with an average intermediate dis- 
tance 10 km. The pressure data were registered at each station with a sample 
time 1 min. An accuracy of measurement was 0.01 %. The recorded signal 
contained obvious diurnal oscillations of the pressure; the average value of these 
oscillations was 0.1 % of atmospheric pressure (~100 Pa), meanwhile, the dif- 
ferential pressure amplitude for neighbour stations had only 10 % accuracy. 

For modeling the angular free-mass interferometer response, we used the 
simple hydrostatic approach calculating the mutual angle between plumb lines 
of suspensions according to the formula 

(22) 

9 v f Fife | 

where pt is the pressure at the z-station, r,fc is the distance between neighbour 
stations (a linear interpolation p(nx) = Kp(x), k < 1 was used). In figure 4 
(a,b), the spectrum of mutual angular fluctuations for pendula separated by the 
distance L is shown. It demonstrates the flicker noise at zero frequencies, diurnal 
peak and a moderate noise level < 10~ n rad-Hz 1/2 at times over few hours. 
In fact, real noise can be less because the noise level in figure 4 is determined 
just by the unsufficient accuracy of the pressure measurements. 
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5 Conclusions 



Our estimations as well as computer simulation of the pendula behaviour in 
real variable pressure field show that, under normal atmospheric conditions, the 
mutual angle between two suspensions spatially separated by 3-4 km fluctu- 
ates with the standard spectral deviation of the order of 1CP 11 rad-Hz -1 / 2 in 
the frequency domain 1CP 3 -j- 1CP 4 Hz, although some rare events (cyclones, 
storms, etc.) can produce one-two order larger perturbations. Such events, 
however, can be forecasted in advance, registered and removed from standard 
data records, or simply the measurements must be eliminated in such times. 
Thus it is quite possible to achieve the desirable accuracy ~ 10~ 13 -j- 1CP 14 rad 
for the measurement time of 10 4 -j- 10 6 sec for detecting the inner core motion. 

Our rough approximations of gravity noises for the turbulent atmosphere 
have shown that the high frequency tail of Newtonian atmospheric fluctuations 
can provide some difficulty for gravitational wave detection at very low frequen- 
cies / < 10 Hz. This point, however, requires a more detailed investigation with 
more adequate models of turbulent currents in the high atmosphere. 
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Figure captions 

Figure 1. Effect of the cyclone dynamics. 

Figure 2. Density d istribution over the cyclone (1), "step" p =const for x < 
xo (la), "circle" p <~ yjl — x 2 Jxq for x < x a (lb), "harmonics" p ~ cos(7tx/.t ) 
for x < x and "Gaussian" p ~ cxp(— x 2 /xq) for all x, angular response as 
function of time (2) and its power spectum (3). 

Figure 3. The pressure spectrum on the Earth surface. Calculated curve 
with h = 18 km, V = 3 m/sec and p = 0.55 kg/m 3 (a) and experimental data 
(b). 

Figure 4. Spectral density of angular response calculated using data of the 
Japan Geographic Institute (Tsukubo Metereological Station). Nord-South di- 
rection (a) and East-West direction (b). 
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